In [4]:
import sys, os
#sys.path.append("C:/Users/lebobcrash/Documents/GitHub/Theano/")
import theano.tensor as T
#sys.path.append("/home/bl3/PycharmProjects/GeMpy/GeMpy")
sys.path.append("C:/Users/lebobcrash/Documents/GitHub/pygeomod/pygeomod")
#sys.path.append("/home/miguel/PycharmProjects/GeMpy/GeMpy")
sys.path.append("C:/Users/lebobcrash/Documents/GitHub/GeMpy/GeMpy")
sys.path.append("/home/bl3/PycharmProjects/GeMpy/GeMpy")
sys.path.append("/home/bl3/PycharmProjects/pygeomod/pygeomod")
sys.path.append("/home/miguel/PycharmProjects/GeMpy/GeMpy")
import GeoMig
import importlib
#import geogrid
importlib.reload(GeoMig)
import importlib
#importlib.reload(GeoMig)
import numpy as np
import pandas as pn
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import matplotlib
import matplotlib.cm as cm
from skimage import measure
os.environ['CUDA_LAUNCH_BLOCKING'] = '1'
np.set_printoptions(precision = 6, linewidth= 130, suppress = True)
from mpl_toolkits.mplot3d import axes3d
import matplotlib.pyplot as plt
from matplotlib import cm
#%matplotlib inline
%matplotlib notebook
In [7]:
GeoMig.Interpolator??
In [ ]:
testmodel = GeoMig.Interpolator(0,140,
0,10,
0,100,
u_grade=0)
# set grid resolution
testmodel.set_resolutions(140,1,100)
testmodel.create_regular_grid_3D()
# compile
testmodel.theano_compilation_3D()
In [79]:
testmodel.load_data_csv("foliations", os.pardir+"/input_data/foliations.csv")
testmodel.load_data_csv("interfaces", os.pardir+"/input_data/interfaces_.csv")
In [80]:
testmodel.formations
Out[80]:
In [81]:
testmodel.set_series({"Initial": (testmodel.formations[0],
testmodel.formations[1],
testmodel.formations[3],
testmodel.formations[4]),
"Unconformity": testmodel.formations[2]
})
In [82]:
testmodel.series
Out[82]:
In [83]:
testmodel
Out[83]:
In [98]:
testmodel.compute_potential_field("Initial", verbose = 1)
testmodel.plot_potential_field_2D(direction = "y", cell_pos = 0, figsize=(7,6), contour_lines = 20)
In [91]:
testmodel.potential_interfaces
Out[91]:
In [92]:
testmodel.compute_potential_field("Unconformity", verbose = 1)
testmodel.plot_potential_field_2D(direction = "y", cell_pos = 0, figsize=(7,6), contour_lines = 20)
In [ ]:
In [99]:
testmodel.compute_block_model??
In [125]:
# Setting the extent
sandstone = GeoMig.Interpolator(696000,747000,6863000,6950000,-20000, 2000,
range_var = np.float32(110000),
u_grade = 9) # Range used in geomodeller
# Setting resolution of the grid
sandstone.set_resolutions(40,40,150)
sandstone.create_regular_grid_3D()
# Compiling
sandstone.theano_compilation_3D()
In [126]:
sandstone.load_data_csv("foliations", "C:/Users/lebobcrash/Documents/GitHub/GeMpy/input_data/a_Foliations.csv")
sandstone.load_data_csv("interfaces", "C:/Users/lebobcrash/Documents/GitHub/GeMpy/input_data/a_Points.csv")
In [127]:
sandstone.set_series({"EarlyGranite_Series":sandstone.formations[-1],
"BIF_Series":(sandstone.formations[0], sandstone.formations[1]),
"SimpleMafic_Series":sandstone.formations[2]})
sandstone.series
Out[127]:
In [128]:
sandstone.compute_potential_field("EarlyGranite_Series", verbose = 1)
sandstone.plot_potential_field_2D(direction = "y", cell_pos = 13, figsize=(7,6), contour_lines = 20)
In [ ]: